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Abstract 

We give a simple solution of the equation {pu')' + qu = uj^u when- 
ever a nontrivial solution of {pu'^)' + quQ = is known. The method 
developed for obtaining this result is based on the theory of pseudoan- 
alytic functions and their relationship with solutions of the stationary 
two-dimensional Schrodinger equation. The final result, that is the 
formula for the general solution of the equation {pu')' + qu = uP'u has 
a simple and easily verifiable form. 

1 Introduction 

The problem of solving the equation 

5a;(p(x)9a;u(a;)) + q{x)u{x) = uj'^u{x) (1) 

by a known nontrivial solution of the equation 

dM^)dxUoix)) + q{x)uo{x) = (2) 

where p, q, u, mq are complex valued functions of the real variable x and u is 
an arbitrary complex constant is of fundamental importance due to numerous 
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situations in mathematical physics where it arises. For example, when the 
method of separation of variables is applied to the equation 

div{PVv) + Qv = 

where P and Q possess some symmetry sufficient for separating variables very 
often one can arrive at the equation ([1]), and it is really desirable to have 
a possibility to solve only one equation ([2]) and to derive from its solution 
the solution of ([1]). Moreover, in many important cases the solution of ([2]) is 
known. For example, consider the conductivity equation 

div{P\/v) = 

and suppose, e.g., that P is a function of one Cartesian variable (for a recent 
work motivating this example see, e.g., [4J). Separation of variables leads to 
the equation 

d^{P{x)d,j;u{x)) = uj'^u{x) 

and the solution of the corresponding equation ([2]) is given, a particular 
nontrivial solution can be chosen as Mq = 1. Thus, to have a method allowing 
us to transform uq into u would mean a complete solution of the original 
problem. 

There are dozens of works dedicated to the construction of zero-energy 
solutions of the Schrodinger equation (see, e.g., [2J, [3j). With the aid of the 
results of the present paper these solutions can be used for obtaining solutions 
for all other values of u. These are just some immediate applications of the 
main result of the present work which can be also applied in spectral and 
scattering theories for the Schrodinger equation as well as for studying the 
Riccati equation. 

Here we give a simple and general solution to the problem of solving ([1]) 
by a solution of ([2]). In obtaining it we used some classical results from pseu- 
doanalytic function theory mainly developed by L. Bers in fourties and fifties 
of the last century as well as some recent results of the author about the re- 
lationship between pseudoanalytic functions and solutions of the stationary 
two-dimensional Schrodinger equation. Although the main result of this pa- 
per was obtained with the aid of the machinery of pseudoanalytic function 
theory it can be easily understood, verified and analized even without those 
tools. The reader not interested in the method derived for obtaining it can 
start reading this paper from equations fH2|) . (H3l) . 
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2 Some definitions and results from pseudo- 
analytic function theory 

2.1 The Vekua equation and pseudoanalytic functions 

The main object of study in the pseudoanalytic function theory is the fol- 
lowing equation 



which is usually called the Vekua equation. Here a, b and W are complex 
valued functions of a complex variable z = x + iy, = djW = ^{dx + idy)W 
and W is a complex conjugate of W. 

Solutions of ([3]) are known as pseudoanalytic (=generalized analytic) func- 
tions. Basic results on them can be found in two well known monographs 
[Ij and [9j. In particular, as a part of pseudoanalytic function theory, L. 
Bers created a well developed theory of formal powers which are solutions 
of ([3]) and locally, near a centre zq, they behave as usual powers {z — zq)'", 
n = 0,1,2,.... Under quite general conditions he (joint with S. Agmon) 
proved the expansion theorem which establishes that any solution of ([3]) in a 
domain of interest Q can be developed into an infinite series of formal powers 
and that this series converges normally (uniformly on any compact subset of 
Q). In posterior works (e.g., [8]) stronger results were obtained guarantee- 
ing the completeness of the infinite system of formal powers in the space of 
pseudoanalytic functions in the sense of the C(fi)-norm. 

In the present work we do not need to consider the general form ([3]) of 
a Vekua equation. Instead, a very special case is sufficient for our purposes. 
In this case the situation with expansion and convergence theorems is even 
simpler. We give corresponding details in subsection 12.31 Here we start with 
some basic definitions and facts from pQ. 

Definition 1 A pair of complex functions F and G possessing in fl partial 
derivatives with respect to the real variables x and y is said to be a gener- 
ating pair if it satisfies the inequality lm{FG) > in Q. The following 
expressions are known as characteristic coefficients of the pair {F, G) 



W^ = aW + bW 



(3) 



a(F,G) 



FG^-FjG 
FG-FG 



hF,G) 



FG^~FjG 
FG-FG 
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FGz — F^G FG z — FzG 

where the subindex 'z or z means the apphcation of dz or dz respectively. 

Every complex function W defined in a subdomain of Q admits the unique 
representation W = (pF + tpG where the functions and ip are real valued. 
Sometimes it is convenient to associate with the function W the function 
w = (j) + iip. The correspondence between W and u is one-to-one. 

For W e C\n) the (F, G)-derivative W = exists and has the 

form 

W = ^zF + i^zG = Wz- A^F,G)W - B^F,G)W (4) 

if and only if 

+ = 0. (5) 
This last equation can be rewritten in the form Q: 

Wz = a(F,G)W + b(F,G)W. (6) 

Solutions of this equation are called (F, G)-pseudoanalytic functions. It is 
said that (F, G) is a generating pair corresponding to the Vekua equation 
(Q. If is (F, G')-pseudoanalytic, the associated function w is called (F, G)- 
pseudoanalytic of second kind. 

Remark 2 The functions F and G are (F, G)-pseudoanalytic, and F = G = 
0. 

Definition 3 Let (F, G) and (Fi, Gi) - be two generating pairs in Q. (Fi, Gi) 
is called successor of{F,G) and (F, G) is called predecessor of{Fi,Gi) if 

a(Fi,Gi) = 0.{F,G) O-nd ^(Fi,Gi) = —B(F,G)- 

The importance of this definition becomes obvious from the following 
statement. 

Theorem 4 Let W be an {F,G)-pseudoanalytic function and let (Fi,Gi) be 
a successor of (F, G). Then W is an (Fi, Gi) -pseudoanalytic function. 
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Definition 5 Let {F, G) be a generating pair. Its adjoint generating pair 
{F, G)* = {F*, G*) is defined by the formulas 

2F 2G 
G* 



FG-FG FG-FG 
The (F, G) -integral is defined as follows 

J Wd(F,G)Z = F{zi) Re J G*Wdz + G{zi) Re J F*Wdz 

where F is a rectifiable curve leading from zq to Zi. 

liW = 4>F + ipG is an (F, G)-pseudoanalytic function where and are 
real valued functions then 



Wd^F,G)Z = W{z) - (t>{z^)F{z) - i^{zo)G{z), (7) 



and as F =G = 0, this integral is path-independent and represents the (F, G)- 
antiderivative of W. 

2.2 Generating sequences and Taylor series in formal 
powers 

Definition 6 A sequence of generating pairs {{Fm, Gm)}, m = 0, ±1, ±2, . . 
is called a generating sequence if {F^^i,Gm+i) is a successor of {Fra,Gm)- 
If {Fq, Gq) — {F, G), we say that {F, G) is embedded in {{F^m Gm)}. 

Theorem 7 Let {F, G) be a generating pair in Q. Let Qi be a bounded 
domain, Qi C Q. Then {F, G) can be embedded in a generating sequence in 

Definition 8 A generating sequence {{Fm, Gm)} is said to have period /i > 
if (Fm+iJ., Gm+n) is equivalent to (F^, Gm) that is their characteristic coeffi- 
cients coincide. 

Let W be an (F, G')-pseudoanalytic function. Using a generating sequence 
in which (F, G) is embedded we can define the higher derivatives of W by 
the recursion formula 

lyioi^ty; m = l,2,.... 

dz 
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Definition 9 The formal power Zm\a, zq; z) with center at zq G fi, coeffi- 
cient a and exponent is defined as the linear combination of the genera- 
tors Fm, Gm with real constant coefficients X, n chosen so that XFmizo) + 
nGm{zo) = a. The formal powers with exponents n = 1,2, . . . are defined by 
the recursion formula 

Zt\a, zo; z) = nf Z^Z;l\a, z,; Qd^FMC- (8) 

This definition imphes the following properties. 

1. Zm\a, Zq] z) is an {Fm, G'm)-pseudoanalytic function of z. 

2. If a' and a" are real constants, then Zm\a' +ia" , zq; z) = a'Zm\l, z^, z) + 
a"zin\i, zo;z). 

3. The formal powers satisfy the difierential relations 

d{F^,G^)Zm\a,Zo;z) _ (n-l), . 
— \(^-,Zo]Z}. 

4. The asymptotic formulas 

Z''m\a,ZQ;z) r-^ a{z - ZqY, z Zq 

hold. 

Assume now that 

oo 

iy(z) = 5^Z(«)(a„,zo;^) (9) 

n=0 

where the absence of the subindex m means that all the formal powers corre- 
spond to the same generating pair {F, G), and the series converges uniformly 
in some neighborhood of Zq. It can be shown [I] that the uniform limit of a se- 
ries of pseudoanalytic functions is pseudoanalytic, and that a uniformly con- 
vergent series of (F, G)-pseudoanalytic functions can be (F, G)-differentiated 
term by term. Hence the function 1^ in ([9]) is (F, G)-pseudoanalytic and its 
rth derivative admits the expansion 

oo 

W^'-\z) = ^n(n - 1) ■ • • (n - r + l)Z("-'-)(a„, zo; z). 
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From this the Taylor formulas for the coefficients are obtained 



(10) 



Definition 10 Let W{z) he a given {F,G)-pseudoanalytic function defined 
for small values of \z — zq\. The series 



(11) 



n=0 



with the coefficients given by OOj) is called the Taylor series of W at Zq, 
formed with formal powers. 

2.3 An important special case 

In the present work in fact we will need the formal powers in the case when 
the generating pair has the form 

iT{y) 



F{x,y) 



aix] 



and G{x,y) 



T{y) a{x) 

where a and r are real-valued functions of their corresponding variables. For 
simplicity we assume that Zq = and F(0) = 1. In this case (see [1]) the 
formal powers are constructed in an elegant manner as follows. First, denote 

x(°)(x) = x(o)(x) = f(°)(2/) = y^'^Hy) = 1 

and for n = 1,2, ...denote 



X(")| 



x] 



for an odd n 



n 



n 



/x("-)(0^ 



X 

y"x("-i)(^)a2(0c^^ for an even 



n 



X("-i)(0fT2(0rf^ for an odd 



n 



X) = i 







n/X("-i)(0^ 





for an even n 



y(n)(^) 



for an odd n 







y 



n J ^\r])T'^{r])d'r] for an even n 



y 

n j Y^"''^\r])T'^(r])dr] for an odd n 



k 

Then for a = a' + ia" we have 



where 



r 



and 



,Z(")(a,0,z) = a 



ia"^Qx("-^)^^F(^) for 



\ %d'Y^{^\x'^''~^\^Y^^ for 



an odd n 



an even n. 



(12) 



(13) 



Consider a rectangular domain containing the origin as an internal point. If 
cr and r are continuously differentiable and bounded together with l/cr and 
1/r on their respective intervals any (F, G)-pseudoanalytic function can be 
represented in the form of a Taylor series in formal powers ffTTl) with = 
and the series converges normally on the domain of interest. 
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2.4 The main Vekua equation 

The equation of the form 

= (14) 

where / is a nonvanishing continuously differentiable real-valued function 
will be called the main Vekua equation. As was shown in [5] this equation is 
related to the stationary Schrodinger equation much in the same way as the 
Cauchy-Riemann system to the Laplace equation. We formulate two results 
which will be used throughout the paper. 

Theorem 11 Let W = Wi+iW2 be a solution of [i4\ )- Then the function 
Wi is a solution of the stationary Schrodinger equation 

- AWi + qiWi = mn (15) 

with qi = Af / /, and W2 is a solution of the associated stationary Schrodinger 
equation 

- AW2 + q2W2 = mn (16) 
where q2 = 2(V/)V/' - Qi and (V/)^ = + /2. 

Notation 12 Consider the equation 

dzV = ^ (17) 

in a whole complex plane or in a convex domain, where $ = $1 + z$2 is a 
given complex valued function such that its real part $1 and imaginary part 
$2 satisfy the equation 

dy^l - 9,$2 = 0, (18) 

then as is well known there exist real-valued solutions (p to equation p7| ) 
which can be reconstructed up to an arbitrary real constant c in the following 
way 

^{x, y) = 2( [ $i(r7, y)dr] + f $2(xo, ^d^] + c (19) 

where {xo,yo) is an arbitrary fixed point in the domain of interest. 
By A we denote the integral operator in fT^) : 

A[^]{x,y) = 2U <^^{r],y)d7]+ T <l>2(xo, e)c?a + c 

\Jxo Jyo J 
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Note that formula / fTPj) can he easily extended to any simply connected domain 
by considering the integral along an arbitrary rectifiable curve T leading from 
{xo,yo) to {x,y) 

(p{x, y) = 2(^j <S>idx + ^2dy^ + c. 

Theorem 13 Let Wi be a real valued solution of ( f73]) in a simply con- 
nected domain VL. Then the real valued function W2, solution of / fl^) such 
that W = Wi + iW2 is a solution of [T^ , is constructed according to the 
formula 

W2 = f-'Aitf'd,if~'W,)). (20) 

Given a solution W2 of l[TS\) . the corresponding solution Wi of l[TS\) such 
that W = Wi + iW2 is a solution of ( fl^ , is constructed as follows 

= -fA{if-^d^{fW2)). (21) 

Remark 14 Observe that the pair of functions 

F = f and G = j (22) 



IS a 



generating pair for [14\ )- 



3 The main result 

Consider the equation 

i-dl + qix))gix)=0. (23) 

We suppose that q and g are real-valued functions and that on some interval 
Ix of the independent variable x there exists a bounded nonvanishing solution 
go e C^iylx) such that 1/go is also bounded. For simplicity we suppose that 
the interval Ix includes the point x = and that goifS) = 1. Our first goal is 
to solve the equation 

{-dl + q{x) ± uj'^)u{x) = 
for any real constant uj. We start with the "+"-case. 
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3.1 The "+"-case 

Consider the equation 

{-dl + q{x) + uj'^)u{x) = 0. (24) 

Let us notice that for the equation 

{-A + q{x)+uj^)U{x,y) = (25) 

where A = + we can immediately propose a particular solution, e.g., 

f{x,y) = go{x)e-y. (26) 

This function does not have zeros on any rectangular domain Q = Ix x ly 
where ly is an arbitrary finite interval of the variable y. For simplicity we 
assume that the origin 2; = is an internal point of the domain Q. According 
to theorem [T3] any solution of (1251) is a real part of a solution of the main 
Vekua equation (1141) where / is defined by (1261) . Moreover, a generating pair 
for this Vekua equation has the form 

Fix, y) = fix, y) = goix)e'^y 

and 

Gix,y) = t/fix,y) = tg^\x)e-^y. 

Now using the results from subsection 12. 31 we can construct the formal powers 
corresponding to equation (fT^ and to this generating pair. We have that 

Z(")(a,0,^) = ^o(a;)e^^Re ,Z'-''\a,0, z) + ig^\x)e-''Hm ,Z(")(a,0,z) 

where ,Z(")(a,0, z) are constructed according to formulas (IT^ and (1131) 
where (t(x) = goix) and r(?/) = e""^^. 

We know that any solution W of f|T^ in Q can be represented in the form 



iy(z) = ^ZW(«n,0,z) 
n=0 

and hence any solution U of (1251) has the form 

00 

?7(x,|/) = ^ReZW(a„,0,;2) (27) 

n=0 

00 

= 5o(2;)e"^5^ (a; Re *Z(")(1, 0, z) + < Re ,Z(")(z, 0, z)) . 



n=0 
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Observe that solutions of are also solutions of (1231) . Consequently, for 
any solution u of there exists such set of real numbers {a^, a^}^o ^^^^ 



u{x) = (7o(a;)e"^5Z 0, 2) + < Re 0, z)) . (28) 

n=0 

In other words there exist such sets of coefficients {aj^, a^}^Q that 

dy ( e"^^ (a;. Re (1, 0, 2) + < Re (i, 0, ;z)) ) = 0. (29) 



n=0 



Obviously, if this condition is fulfilled, the resulting function (1251) is a solution 
of dH. 

Let us analyse equation (|29|) . First of all we have that for an odd n, 



Z(")(z,0,z) = 



i=o 

and for an even n, 



i=o 

Thus, we obtain for an odd n, 



Z(")(l,0,z) = ^ 



Re,Z(")(l,0,z)= ('^^("-^Vf^^) 



even j=0 



Re.Z(")(z,0,z)= 

odd j=l ^ 
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and for an even n, 



Re*z(")(i,o,z) = J2 ('!)^^""'Vy^^), 

even j=0 ^ 



Re,Z(")(z,0,^)= Yl C")^ 

odd 7=1 ^-^^ 



It is somewhat more convenient for what follows to rewrite these formulas in 
the following equivalent form 



Re *z(")(i,o,^) = 

odd k=l 

n 



for an odd n (30) 



even k=0 



and 



Re *ZH(l,0,z) = Y 

even fc=0 
n 

Re *Z(")(i,0,^) = Y 

odd fe=l 

We remind that 

f y 

I J F("-^)(77)e2'^^d77 for an odd 



y 

I j Y^'^~^\r])e~'^'^^dr] for an even 



for an even n (31) 



yW(y) = < 



n 



n 



yW(y) = < 



n 



n 



y(n i)(^7^)e '^'^'^d'q for an odd n 
y («■-!) (,y)e2'^^d77 for an even n 
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Thus, from (1501) and ([21]) we have 

n 

dy Re 0, z) = ^ (n - A;) Q)x('=)i"-'=F("-'=-i) 

"'^^ for an odd n 



even k=0 

and 



(32) 



aj^Re ,Z(")(l,0,z) = J] (n - A;) (^)X('=)i"-'=F("-'^-i) 

even fc=o f^j, even n > 



odd fc=l 



(33) 

Now returning to equation fl29|) we observe that it is equivalent to the equa- 
tion 

oo 

J2K Re (1, 0, z) + dy Re (1, 0, z)) 

n=0 

+< (cu Re (^, 0, z) + dy Re (^, 0, 2)) ) = 

from which using (1301) . (!3T|) and (!32|) . (!33|) we obtain that (l29l) can be written 
as follows 

even n=2 even A;=0 

+< ^ ('^')xW(^?^""'^ + - A;)y("-'=-^)e-2"^)) 

odd k=l ^ ^ 

odd n=l odd k=l ^ ' 

+ < ^ 2"^'=+^ r^^J (u;?^"-'^) + (n - A;)F(""'=-i)e~2"^)) = 0. (34) 

even fc=0 

In order that this equality hold identically the expressions corresponding to 
different X^^'> and for all n should vanish identically. Combining all 
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terms multiplied by X^^^ we obtain the equation 

oo 

even n=2 

oo 

+ Yl +^?^""'^e-2"^) = 0. (35) 

odd n=l 

Gathering all terms multiplied by X^^^ we obtain the second equation 

oo 

J2 <^"n(a;F("-^) + (n - i)y("-2)e-2-2/) 

even n=2 
oo 

odd n=l 

which can be rewritten as follows 

oo 

even n=2 

oo 

+ J2 + + nf ("-i)e~2-?') = 0. (36) 

odd n=l 

Gathering all terms multiplied by X^'^\ X^^\ ... we obtain an infinite system 
of equations which fortunately we do not need to solve. Here we are reasoning 
along the following lines. First of all we observe that if such sets of coefficients 
{a^, a"}^o 6xist that all the equations derived from flM|) are satisfied, they 
do not depend on functions X^^^ and X^'^^ but only on y^"-* and Y^^\ Thus, 
they can be constructed independently of the concrete form of go and hence 
of the potential q. Second, we know that such sets of coefficients exist. This 
is due to our earlier observation that solutions of fl2^ are also solutions of 
fl25|) and hence they can be written in the form fl27|l . 

These two arguments lead to the following surprising solution. We can 
take any q, for example, q = and any pair of independent solutions of the 
resulting Schrodinger equation ( 12^ and to obtain their corresponding sets 
of coefficients. These two sets will be universal in the sense that the general 
solution of fl2^ for any other q will be constructed with the aid of this pair 
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of sets of coefficients just changing the generating function and obtaining 
a corresponding system of functions X^""^ and X^"'\ On the first glance this 
conclusion can appear against the intuition, nevertheless its more detailed 
analysis as well as the final result convince that it is really natural. Thus, in 
the next subsection we construct such pair of sets of Taylor coefficients (in 
formal powers). 

3.2 Two sets of Taylor coefficients 

Here we consider the case g = 0. Then the Schrodinger equation fl2^ becomes 

i-dl + uj'Mx) = 0. (37) 

Note that the corresponding equation (!23|) has the form d'^g{x) = and 
possesses a suitable particular solution satisfying all the requirements (see 
the beginning of section E]) (^o = 1- Then / = e'^^ and the main Vekua 
equation in this case has the form 

= y TF. (38) 

Let us take two independent solutions of (!37|l u~^{x) = e'^^ and u^{x) = e~'^^. 
First, we obtain the set of coefficients {a^,a^}^o function The 

first step consists in constructing the corresponding conjugate metaharmonic 
function (see theorem [T3|) 

v+ = e-'^^Z(ie2'^^5^e'^(^-^)) = e-'^M(-(l + i)e^(^+^)). 

2 

We have 

A(|(l + i)e^(^+^)) = e^(^+^) + c. 

We choose c = 0, then v"^ = e^^. Thus, one of the solutions of the main 
Vekua equation (!38|) such that = e^^ is its real part has the form 

W+ = {1 + (39) 

Now, in order to construct its corresponding Taylor coefficients (in formal 
powers) we notice that Ai^p^Q-^ = and B(^pQ) = —iuj/2 (here F = e'^^ and 
G = ie~'^y). Thus, the operation of the {F, G)-derivative has the form 

iuj — 

w = w, + —w. 
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For the function (1391) we have 



and it is easy to see that the ra-th (F, G')-derivative of has the form 

We obtain that the Taylor coefficients (in formal powers) of the function (13^ 
at the origin have the following simple form 

.71 

«;l = ^(i + 0- (40) 

nl 

In a similar way we study the case of the function . The correspond- 
ing pseudoanalytic function W~ has the form W~ = (1 — i)e~'^^, and the 
corresponding Taylor coefficients at the origin are as follows 



'-cur 



l-i). (41) 



Let us notice that from the fulfillment of ( l35i) with the coefficients of the 
form (Ho!) or (HTj) there follows the fulfillment of (136!) and of all subsequent 
equations corresponding to X'^'^\ X^'^\ etc. This is because of the fact that 



3.3 General solution of ( l24| ) 



Now with the aid of the sets of coefficients pUj) and (HTj) we proceed in 
obtaining the general solution of (!24|) with any potential g for which a solution 
(70 of ( j23l) satisfying the nonzero and boundedness requirements exists. From 
( I28l) we have that the general solution of ( |M1) has the form 



U = CiUi + C2U2 

where ci and C2 are arbitrary real constants and ui, U2 are defined as follows 
ui{x) = (7o(a;)e"^y^ (Re ,^("^(1, 0, 2) + Re ,Z^''\i,0, z)) 

n=0 
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and 



CO , N ri 

U2{x) = (7o(a;)e"^ V^^^^ (Re 0, z) - Re 0, z)) 

n=0 

which according to ( l30l) and ( l3Ti) can be written in the following form 

E ^( E .'-'■(nx<"i'<"-" 

/en n=0 even fc=0 
_^ ^ ^n-fc+l/''^^J^(fc)y{n-fc)^) 



even n=0 even fc=0 



odd fe=l 
oo 



+ E ^( E i"--(^)x(-)F<'.-') 

odd ra=l ■ odd k=l ^ ^ 



even A;=0 



and 



«2(a;) = ^7o(x)e-n E ^( E ^^'M 

even n=0 even fc=0 

odd k=l 

- E ^( E r-'(;;)x«y'-'> 

odd n=l odd A:=l ^ ^ 



even fc=0 



As we know that both expressions are independent of y in order to simplify 
them we can substitute any value of y. Of course, the easiest way is to 
substitute y = because by definition all F(")(0) and y(")(0) for n > 1 are 
equal to zero, and y(°)(0) = F(°)(0) = 1. Thus, finally we obtain 



n! — ' w. 

even n=0 odd n=l 
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and 



where 



even n=0 



odd n=l 



X = < 



X(")| 



x = < 



n y X("-i)(0^o(0^^^ for an odd n 



X 

n j X^''^~^\^)gQ'^{^)d^ for an even n 



n j X("-i)(0^o"^(Of^^ for an odd n 



n j X("-i)(0^o(0^^^ for an even 



(43) 

(44) 



(45) 



(46) 



n 



In the next subsection we vahdate this resuh by a direct substitution into 
equation 



3.4 Validating the result 

In order to substitute and or equivalently 



v,{x)=go{x) J2 



even ri=0 



and 



odd n=l 

into equation we first make some helpful observations. 

It is well known that a nonvanishing solution go of (1231) allows us to 
factorize the Schrodinger operator as follows 



^0 



^0 



(47) 
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The first order operators in their turn can be factorized as well, so we obtain 



Now let us consider vi. By definition, for an even n we have 

X 







alar 



Thus, application of the operator d1 — q to QqX^'^^ for an even n and n > 2 
(for n = the result is zero) gives us 

{dl - q) (^?oX(")) = g^^d^gld^X^^^ = ng^'d^X^^^'^ 
= (n-l)n(7oX("-2). 

Then 

even n=2 ^ ^ 

oo n 

even n=0 

In a similar way one can verify that V2 is a solution of (!24|) as well. Note 
that according to the general result formulated in subsection 12.31 both series 
in Vi and V2 are uniformly convergent on the interval Ix- This fact can be 
quite easily verified as well by estimating the integrals in X*^"^ and X*^") by the 
supremum of the functions g^ and g^^ multiplied by successive antiderivatives 
of X. 

3.5 The "-" case 

Consider the equation 

{-dl + q{x) - uj'^)u{x) = (48) 
and the corresponding two-dimensional equation 

(-A + g(a;) -cj2)t/(x,y) = 0. (49) 
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Its particular solution can be chosen as 

fix,y) = goix) cos ujy 

which is different from zero on the rectangular domain Q = x (~^, 
In order to obtain the general solution of ( HHl) in fact we should only obtain 
two sets of Taylor coefficients as in subsection 13.21 For this, once more we 
take g = and consider two linearly independent solutions of the equation 
{dl + u'^)u{x) = 0, u'^{x) = cosujx and u^{x) = sinux. The next step is to 
construct and v~. We have 

, 1 — / 2 o /cosciJx\\ 

= A I cos ujyOz = — sm ujxianuy 

cos ujy \ \ cos ujy J J 

(we have fixed the arbitrary constant as zero). Thus, 

= cos ujx — i sin ujx tan uy. 

In a similar way we obtain 

W~ = sin UJX + i cos ux tan uy. 

Noting that the definition of the {F, G)-derivative in this case has the form 

iui — 
W = W^- —tanujyW 

we obtain the following relations = —ujW~ and W ~ = ujW~^ and hence 
the following formulas for the corresponding Taylor coefficients in formal 
powers in the origin 

~ ^ even n and = for an odd n, 

nl 

a~ = for an even n and a~ = ; — for an odd n. 

n! 

Thus we arrive at the following general solution of equation ( 148 p for any 
potential q admitting a particular solution with the described above prop- 
erties 

U = CiUi + C2M2 
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with 

oo 

[lUJ 



nl 

even n=0 



oo ■/ ■ \ 



and 

nl 

odd n=l 

where X(") and are defined bv (l44l)-fH6|). 
3.6 Complex potential 

It is clear that the results obtained in the preceding subsections remain valid 
in the case of a complex valued potential q and an arbitrary complex number 
u. Consider the equation 

{-dl + q{x) + uj^)u{x) = (50) 

where q and u are complex valued and u is any complex number. We assume 
that go is a nonvanishing solution of the equation {—dl+q{x))go = satisfying 
the boundedness requirements. Then the general solution of fl50l) has the form 

U = CiMi + C2U2 (51) 

where Ci and C2 are arbitrary complex constants and Mi, U2 are defined as 
follows 

00 „ 

^^ = 90 E ^^XW (52) 

even n=0 

and 

00 „ 

«2 = ^70 E (53) 

odd n=l 

where as before X*^") and X*^") are defined by 



3.7 Equation {pu')' + qu = u^u 

Having obtained the general solution of fISUl) it is easy to obtain the general 
solution of the more general equation 

dxipdxu) + qu = uP'u (54) 
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where p G C'^{Ix) is a nonvanishing complex valued function, q and u satisfy 
conditions from the preceding subsection. We assume that 

dx{pdxgo) + qg^ = Q (55) 
and observe that the following factorization holds 



where (7 = p^^'^go and by analogy with fl5Tl) - fl53l) we obtain the following 
solution of 

U = CiUi + C2M2 

where 



wi = go 



^ n\ 



Hn) 



even n=0 



and 



odd n=l 

where the definition of X^"'^ and is slightly modified: 



(56) 
(57) 



(58) 



(59) 



X^^'^KOgiiOd^ for an odd 



n 



x) = < 







n I X*^" "'^HOs' ^(0*^^ even n 









X 

j X^''-^\i)g^^{i)di for an odd 


X 

j X^''-^\i)gl{Cjdi for an even 



n 



n 



(60) 



(61) 
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Let us verify that Ui is indeed a solution of We have 



{dxpd^ + q)ui = go ^d^ g'^d^ ^ —X 



"(n) 



even n=0 



oo „ 
UJ 



0-1 V — - 9x("-^) 

(n-2)!" 



n - 1 ! 

even n=2 



even n=2 

In a similar way the solution U2 can be verified as well. 

4 A remark on spectral problems 

First, let us notice that at least in the case of real- valued coefficients p and q 
under the conditions that p, p' and q are continuous an appropriate nonvan- 
ishing particular solution go of (!55|) always exists. In this case the equation 
{pg'Y + qg = possesses two linearly independent solutions gi and g2 whose 
zeros do not coincide. Then g^ can be chosen as follows go = gi + ig2- Thus, 
the proposed solution can be obtained in a quite general situation 

and as we show in this section the corresponding spectral problems reduce 
to the problem of finding zeros of related analytic functions. 

Sometimes it is slightly more convenient to consider the function fl58l) 
divided by uj, that is 

°° n— 1 

nl 

odd n=l 

Then for the solutions of fl54|) ui and U2 defined by (|57|) and fl62|) we obtain 
the following equalities 

"i(0) = go(O), ii',(0)=gi(0), (63) 

Now consider a spectral problem for (1541) on the interval Ix = (0,1). For 
example, 

n(0) = and u(l) = 0. (65) 
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Due to the first boundary condition the constant Ci in fl56l) should be chosen 
as zero. Then the spectral problem reduces to finding such values of uo that 

oo 

■U2(l) = fi'o(l) ^^^^r"^*'"''(l) vanish. In other words, this spectral problem 

odd n=l 

reduces to the calculation of zeros of the complex analytic function 



oo 
m=0 



where defined as follows 



for an odd m 

^i^^M — -iM for an even m. 

Note that in some cases the function n will be even entire. For example, when 
the considered Sturm-Liouville problem is regular it is well known (see, e.g., 
[6]) that there exists an unboundedly increasing sequence of its eigenvalues 
which means that there exists an unboundedly increasing sequence of zeros 
of the function /t. 

Let a and (3 be arbitrary real numbers. Consider the following more 
general boundary conditions 

m(0) cosa + m'(0) sina = (66) 

M(a) cos/3 + M'(a) sin/3 = (67) 
together with equation (j54l) . Taking the solutions ui and U2 defined by ( !57ll 



and (1621) and using (]63i) . (l64l) we obtain from mm the following equation for 

ci and C2 

sin a 

ci (5-0(0) cosa + 5(0(0) sin a) + ^'^'^J^y^ = ^ 
which gives the relation 

C2 = 7C1 for a Tin 
where 7 = -50(0)^(0) (50(0) cota + 50(0)) and 

Ci = for a = Tcn. 
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In this last case we arrive at a similar result as in the example considered 
above, thus let us consider the case a ^ Tin. From the definition of Ui and 
U2 we have 



1 = + — 



even n=2 



{n-l)\ 



and 



h = —^2 + > -r^^ 

y uneven n=0 

Then the boundary condition ( 1671) implies the following equation 
(^7o(a)cos/3 + Sin/5) ^ i^x(")(a)+7 ^^X^")' 

sin /? 



v,even n=0 



odd n=l 



go{a)p{a) 



^ (ra-l)! ^ n! / 



^^even n=2 ^ ' even n=0 

Thus the spectral problem flMl) . fl^ . (jnZD reduces to the problem of finding 
zeros of the analytic function 



= E 



m=0 



where dm cLFG defined as follows 



«o = (fl'o(a.) cos/5 + g'Q{a) sin/5) (1 + 7X^^^(a)) + 



7 sin /5 
9o{a)p{a) 



and 



r 



(c/o(a) cos/5 + ^o(a) sm/?) ( + 7 (^+1)1 



sinff / X('"-^)(a) _|_ X('")(a) 
go{a)p(a) \ (m— 1)! ' m! 



for an odd m 



for an even m > 0. 
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5 A remark on the Darboux transformation 

The Darboux transformation is a very useful and important tool studied in 
dozens of works (see, e.g., [7]). It is closely related to the factorization of the 
Schrodinger operator (H7j) . Consider the equation 



9oJ V 9oJ 



Applying the operator — ^ j to both sides and denoting v = (^d^ — ^ j u 
one obtains that solutions of equation fISUl) are transformed into solutions of 



another Schrodinger equation 



which can be written also as follows 



{-dl + r{x) + v{x) = 0, (68) 

where r = 2 j — q. Now, as we are able to construct the general solution 

of (jSOl) by a known solution of (!23l) we can also obtain an explicit form of 
the result of the Darboux transformation. First, let us apply the operator 

(^dx — = godxdQ^ to Ui defined by flS21) . We have 



I \ °° n 

On / ^-^ nl 

' even n=0 

oo „ oo 



^° ^ (n- 1 ! go ^ n\ 

evenn=2^ ' ^^oddn=l 



and in a similar way we obtain 



Qo J Qo — nl 

/ ^^^^ 



Thus, the general solution of the Schrodinger equation fl68l) obtained from 
(|50|1 by the Darboux transformation has the form 



oo „ oo 



even n=0 odd n=l 

where X'^^^ and X*^"^ are defined by 
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6 Conclusions 



The main result of this work allows us to find the general solution of ([T]) 
by a known solution of ([2]) as follows. By formulas (l59l) - (!6T!) the functions 
X^""^ and X^"^ should be constructed and then the general solution of ([T]) 
has the form flS^ - fl551) . An important feature of these formulas consists in 
the fact that the form of X*^"^ and X*^") does not depend on cu. For given 
coefficients p and q they should be calculated only once and then the general 
solution for any uj is obtained by multiplying them by corresponding powers 
of uj. This is very convenient for numerical calculations. Moreover, our 
numerical experiments show that X*^") and X(") can be calculated up to high 
indices with a remarkable accuracy, thus we expect that among other possible 
applications the results of this work will be useful in numerical solving of a 
wide class of boundary value and spectral problems of mathematical physics. 
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